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ABSTRACT 

Elastic-viscoplastic  constitutive  equations  based  on  a single 
internal  state  variable  which  is  a function  of  plastic  work  are 
used  to  calculate  the  response  of  copper  to  a six  decade  change  of 
strain  rates  over  a range  of  temperatures.  Calculations  were  per- 
formed for  the  conditions  of  an  experimental  program  on  copper  by 

Senseny,  Duffy,  and  Hawley,  namely,  temperatures  ranging  from  77 
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to  523°K  and  strain  rate  jumps  from  2 x 10  sec  to  3 x 10  sec 
at  three  strain  levels.  The  computed  results  are  in  good  agreement 
with  the  experiments  and  show  similar  strain  rate  and  strain  rate 
history  effects.  Relations  are  obtained  for  the  temperature  depen- 
dence of  certain  parameters  in  the  equations  which  indicate  corres- 
pondence between  plastic  working  and  temperature  and  between  strain 


rate  sensitivity  and  temperature. 
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INTRODUCTION 

Recent  formulations  of  elastic-viscoplastic  constitutive 
equations,  e.g.  [1,2, 3,4]  propose  to  represent  rate  dependent 
material  response  for  a wide  range  of  loading  histories.  A com- 
mon approach  in  most  of  these  theories  is  that  all  aspects  of 
non-reversible  deformation  are  lumped  into  a single  equation 
relating  the  inelastic  deformation  rate  to  stress,  temperature, 
and  one  or  more  internal  state  variables.  When  this  is  taken 
together  with  the  elastic  deformation  rate  equation,  the  total 
deformation  rate  can  be  expressed  as  a function  of  state  vari- 
ables and  the  current  geometry. 

Most  of  the  recent  work  on  viscoplastic  constitutive  equa- 
tions has  been  motivated  by  the  need  to  develop  suitable  creep 
equations  under  conditions  of  fluctuating  load,  strain,  or  tem- 
perature. The  form  of  the  equations  and  the  manner  of  including 
temperature  dependence  are  therefore  similar  to  high  temperature 
creep  equations.  Other  formulations,  in  particular  that  of  Bodner 
and  Partom  [1,5,6],  were  initially  directed  towards  representing 
strain  rate  effects  over  a wide  range  of  conditions.  Those  equa- 
tions serve  to  calculate  stress-strain  curves  for  imposed  constant 
and  varying  strain  rates.  Although  not  specifically  intended  for 
creep,  the  equations  can  also  represent  creep  and  stress  relaxa- 
tion behavior  [1]. 

The  intention  of  the  investigation  described  in  this  paper 


was  to  examine  the  capability  of  the  constitutive  equations  of 
Bodner  and  Partom  to  represent  a relatively  rate  insensitive  metal, 
copper,  over  a wide  range  of  strain  rates  and  temperatures.  These 
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equations  are  based  on  a single  parameter  to  characterize  the  micro- 
structural  state  with  respect  to  inelastic  deformation.  A recent 
experimental  program  by  Senseny,  Duffy,  and  Hawley  [7]  provides 

shear  stress-strain  curves  for  annealed  (640°K)  OFHC  copper  under 

-4  -1  2 -1 

static  (2  X 10  sec  ) and  dynamic  (3  x 10  sec  ) conditions  for 

temperatures  ranging  from  77°K  to  523°K.  Tests  were  also  conducted 
in  [7]  for  rapid  changes  from  the  low  to  the  high  strain  rate. 
Similar  experiments  have  also  been  conducted  by  Eleiche  and  Camp- 
bell [8]  using  differently  prepared  copper  specimens.  A number 
of  other  investigators  have  obtained  static  and  dynamic  stress- 
strain  curves  for  OFHC  copper  with  comparable  results.  For  con- 
sistency, it  was  decided  to  use  a single  set  of  results,  namely 
[7] , as  the  reference  for  this  study. 

The  steady  strain  rate  curves  of  [7]  were  used  as  the  basis 
for  determining  the  material  constants  in  the  constitutive  equations 
at  the  various  temperatures.  With  these  constants,  response  curves 
for  rapid  changes  of  strain  rate  were  calculated.  Their  relation- 
ship to  corresponding  test  results  was  taken  as  a measure  of  the 
applicability  of  the  equations.  The  equations  have  also  been  used 
to  predict  creep  response  under  constant  and  varying  load  condi- 
tions and  this  will  be  reported  in  a subsequent  paper.  For  this 
case,  a recovery  term  was  added  to  the  basic  equations. 

The  empirical  determination  of  the  material  constants  led  to 
relatively  simple  analytical  forms  for  their  temperature  variation. 
These  should  be  useful  for  calculations  involving  varying  tempera- 
ture conditions  and  in  identifying  the  operative  microscopic  de- 
formation mechanisms. 
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DETERMINATION  OF  MATERIAL  CONSTANTS 


Copper  is  relatively  insensitive  to  rate  of  strain  effects  and 
exhibits  fairly  strong  strain  hardening.  The  experiments  reported 
in  [7]  cover  six  decades  of  strain  rate  and  show  strain  rate  and 
strain  rate  history  effects  which  are  emphasized  in  the  tests  in- 
volving rapid  changes  of  strain  rates.  Rather  than  hypothesizing 
the  manner  of  including  temperature  dependence  in  the  constitutive 
equations,  an  empirical  approach  was  followed.  That  is,  numerical 
exercises  were  undertaken  to  find  suitable  values  for  the  material 
constants  so  that  the  equations  would  represent  the  experimental 
results.  An  objective  in  these  exercises  was  to  minimize  the 
number  of  constants  that  were  temperature  dependent. 

The  basic  equations  in  the  formulation  of  Bodner  and  Partom 
[1]  relate  the  reversible  and  non-reversible  components  of  the 
deformation  rate,  which  are  taken  to  be  always  non-zero,  to  state 
quantities.  For  small  strains,  the  reversible,  i.e.  elastic,  de- 
formation rate  is  given  by  the  time  derivative  of  Hooke's  Law. 

Since  the  tests  used  for  reference  were  in  pure  shear,  only  the 
shear  modulus  G is  required.  The  temperature  variation  of  G is 
small  over  the  range  of  the  tests  and  was  not  significant  in  the 
numerical  calculations.  At  room  temperature,  G for  copper  is 
45000  MPa. 


A general  form  for  the  non-reversible,  i.e, 
formation  rate, which  is  equal  to  the  strain  rate 


inelastic 
for  small 


de- 

strains , 


(1) 
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where  are  the  elastic  stress  components,  T is  the  temperature, 

and  Zj^  are  one  or  more  internal  state  variables.  From  the  flow 
law  of  classical  isotropic  plasticity. 


d? . = As.. 
11  11 


it  follows,  upon  squaring,  that 


Df  = A"J2  (3) 

where  is  the  second  invariant  of  the  plastic  deformation  rate 
and  J2  is  the  second  invariant  of  the  stress  deviator,  A funda- 
mental postulate  in  the  formulation  of  [1,5,6]  is  that  plastic 
deformation  is  governed  by  a continuous  relation  between  and 
J2  involving  T and  Zj^.  Solving  for  A in  (3)  and  substituting  into 
(2)  then  gives  an  equation  for  d^^  in  the  form  of  (1) 

d?.  = [Fi  (J2  ,T,Z,  )/J2  ] '/"s.  (4) 

11  ^ 11 

The  expression  chosen  for  Fi  was  motivated  by  equations  pro- 
posed to  relate  dislocation  velocity  to  stress,  namely, 

of  = F,(J2,T,Zj^)  = Doexp[-(jZ^)‘^(!^)/J2]  (5) 

where  some  of  the  constants  are  temperature  dependent.  The  para- 
meter Z in  Eq.  (5)  can  be  interpreted  as  an  internal  state  variable 
which  provides  a measure  of  the  overall  resistance  to  plastic  flow. 
It  has  the  dimension  of  stress  and  can  be  referred  to  as  the 
"hardness".  As  such,  it  is  similar  to  the  hardness  state  variable 


of  Hart  [2]  but  is  treated  in  a different  manner.  Eq . (5)  is  shown 

in  Fig.  1 for  various  Z values. 
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In  Eq.  (5) , Do  is  a scale  factor  and  n relates  to  strain  rate 
sensitivity  with  higher  n meaning  lower  sensitivity.  Additional 
constants  arise  in  the  evolution  equation  for  Z which  was  taken 
to  be  a function  of  plastic  work  VJ  . In  representing  titanium 

Ir 

for  various  uniaxial  strainina  histories  [1] , the  function 


Z = Zj  -r 


( Z 0 -Z ; ) exp ( -mW 


(6) 


which  introduces  three  more  constants,  was  adequate.  This  was  not 
the  case  for  copper  which  exhibits  strong  strain  hardening, so  m 


was  made  a function  of  W 


P' 


m = mo 


+ m^expi-aW 


(7) 


i' 


which  adds  two  more  constants.  It  is  likely  that  fewer  than  five 
constants  are  actually  needed  for  Z (W^)  for  copper,  Fig.  2,  by 
other  choices  of  the  basic  functional  form  (6) , e.g.  a polynomial, 
but  this  was  not  attempted  here. 

Consideration  of  temperature  dependence  in  the  equation  for 


i 


the  inelastic  strain  rate  has  been  studied  extensively  under  creep 
conditions  and  suggestions  on  similar  relationships  have  been  pro- 
posed for  steady  plastic  flow  at  moderate  temperatures.  High  tem- 
perature (T  > creep  equations  are  generally  presented  in  the 

form  of  (1)  in  which  the  temperature  effect  appears  as  a simple 
multiplying  factor  (1/kT)  where  k is  Boltzmann's  constant.  A gener- 
al form  for  temperature-strain  rate  correspondence  was  proposed  in 
[9],  namely  that  steady  rate  stress-strain  curves  depended  on  the 
single  parameter  eexp(flH/RT).  This  correspondence  was  shown  in 
[9]  to  be  applicable  for  pure  aluminum  over  a range  of  constant 
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strain  rates  and  temperatures.  More  recently,  Campbell  [10]  pre- 
sented a detailed  discussion  on  temperature-strain  rate  relations 
in  metal  plasticity.  He  developed  expressions  for  the  inclusion  of 
temperature  effects  in  strain  rate-stress  equations  based  on  models 

t 

for  deformation  mechanisms.  Some  of  these  expressions  could  be 
expected  to  apply  for  steady  rate  conditions  but,  as  Campbell  also 
emphasizes,  may  not  be  applicable  for  histories  involving  changing 
strain  rates. 

Following  the  extensive  work  on  creep  and  the  correspondence 
law  proposed  in  [9],  it  seemed  reasonable  to  consider  Do  and  n in 
(5)  to  be  functions  of  temperature.  Initially  a number  of  numerical 
exercises  were  performed  on  this  basis,  but  the  resulting  functions 
for  Do (T)  and  n(T)  were  not  reasonable,  e.g.  Do  varied  with  tempera- 
ture by  a few  orders  of  magnitude.  Further  exercises  based  on 
varying  only  Z(Wp)  and  n with  temperature  led  to  the  conclusion 
that  the  function  Z (W^)  could  be  considered  generally  applicable 
with  the  initial  state  Z^  being  temperature  dependent.  Applicability 
of  a single  Z(W  ) curve  implies  a temperature-hardness  correspondence 
in  the  constitutive  equations.  Lowering  the  operating  temperature 
would  be  equivalent  to  increasing  the  initial  hardness  value  which 
could  have  been  achieved  by  prior  work  hardening.  One  of  the  re- 
sults reported  in  [9]  does  imply  this  effect.  Since  is  the 
stable  hardness  value  at  a given  temperature,  the  recovery  (anneal- 
ing) equation  of  hardening  due  to  plastic  work  would  have  the  form 


- Z 


i 
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F2 [z (Wp) 


(8) 
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obtaining  the  material  constants  by  matching  calculated  and 
experimental  stress-strain  curves  could  be  done  on  a trial  and 
error  basis  or  more  formally  by  computations  using  a least  squares 
procedure.  In  general,  a pair  of  stress-strain  curves  at  two  dif- 
ferent steady  rates  are  necessary  and  sufficient  to  provide 
essentially  unique  values  for  all  the  material  constants  at  a 
given  temperature.  Some  variation  in  the  values  would  result  from 
different  criteria  of  curve  matching.  Accuracy  of  matching  curves 
was  attempted  to  within  the  experimental  accuracy  which  was  given 
as  5%  in  [8]  and  was  probably  similar  in  [7]. 

Since  the  experimental  values  of  the  initial  stress  increment 
due  to  the  ]ump  in  strain  rate,  At  (Fig.  1),  should  be  more  exact 
than  the  gap  between  the  two  steady  rate  curves,  it  would  seem  pre- 
ferable to  use  the  (At) 's  in  conjunction  with  the  lower  rate  curve 
to  obtain  the  material  constants.  However,  the  procedure  of  forcing 
agreement  of  the  respective  At  values  tends  to  lead  to  overall 
distortion  of  the  steady  rate  curves.  A number  of  numerical  exer- 
cises were  performed  using  this  approach  but  were  discontinued. 

The  material  constants  obtained  at  the  various  test  tempera- 
tures are  listed  in  Table  1.  As  mentioned,  these  are  not  precisely 
unique  values  and  limited  variations  of  combinations  of  the  con-^tants 

lead  to  almost  identical  results.  The  function  Z (W  ) based  on  these 

P 

constants  is  shown  in  Fig.  2 with  the  values  of  at  the  test  tem- 
peratures indicated.  Calculated  steady  rate  stress-strain  curves 
are  shown  in  Figs.  3-6  with  corresponding  test  results  from  [7). 
Superimposed  on  these  curves  are  the  respective  results  for  rapid 
changes  of  strain  rate.  A detailed  "blow  up"  of  one  of  the  calcu- 
lated jump  curves  is  shown  in  Fig.  7. 
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As  noted  in  Table  1,  the  material  constants  that  were  found 
to  vary  with  temperature  were  and  n.  Plotting  Z^  against  the 
reciprocal  of  absolute  temperature  leads  to  a straight  line.  Fig. 

8,  which  can  be  expressed  analytically  as 

Z^(riPa)  = (5000  + 15T)/T(°K)  (9) 

A plot  of  n directly  against  T,  Fig.  9,  suggests  a linear  relation- 
ship which  can  be  expressed  as 

n = - 0.025T(°K)  + 16.6  (10a) 

Although  Eq.  (10a)  provides  a good  fit  to  the  derived  n values,  it 
is  physically  unrealistic  outside  the  range  of  test  temperatures 
since  it  is  expected  that  n-*-«>  as  T-*-0  and  n should  be  small  but  non- 
zero as  A hyperbolic  or  piecewise  linear  relationship  for 

n(T)  would  be  more  meaningful.  A least  squares  linear  fit  of 
n(l/T)  gives 

n = (835  + 5.0T)/T(°K)  (10b) 

with  fairly  large  deviations.  These  deviations  appear.  Fig.  9,  to 
be  greater  than  expected  from  experimental  scatter  and  variations 
in  curve  matching  criteria,  but  a hyperbolic  relation  such  as  (10b) 
may  be  a reasonable  approximation  over  a large  range  of  T. 

Recovery  of  hardness  (annealing)  was  not  considered  in  calculat- 


ing the  stress-strain  curves  but  would  be  important  under  creep  con- 
ditions. In  fact,  these  constitutive  equations  with  a recovery  term 
added  would  show  the  generally  accepted  result  that  secondary  creep 


is  a balanced  condition  between  hardness  increases  due  to  strain- 
ing and  decreases  due  to  recovery.  An  equation  for  hardness 
recovery  suggested  in  the  literature  in  the  form  of  (8)  is 

Z = - C(Z-Z^)^  (11) 

with  r an  integer.  For  a particular  r,  C could  be  specified  so  that 
the  equations  would  lead  to  a prescribed  constant  creep  rate  at  a 
particular  stress  and  temperature. 

Inclusion  of  the  recovery  term  (11)  also  enables  calculation 
of  transient  effects  during  creep  due  to  increments  or  decrements 
of  load  or  temperature.  An  interesting  phenomenon  is  the  delay 
time  consequent  to  a small  load  decrement,  i.e.,  a period  of  essen- 
tially zero  strain  followed  by  creep  at  the  lower  stress  level. 
Preliminary  calculations  indicate  that  the  equations  do  show  this 


effect. 
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DISCUSSION  OF  RESULTS 

A comparison  of  the  calculated  and  experimental  results , 

Figs.  3-6  and  Table  2 (with  symbols  defined  in  Fig.  1) , leads  to 
the  tentative  conclusion  that,  at  least  for  copper,  viscoplastic 
constitutive  equations  based  on  a single  internal  state  variable 
which  is  a function  of  plastic  work  are  adequate  to  represent  strain 
rate  history  effects  associated  with  rapid  changes  of  strain  rate. 

In  Figs.  3-6  and  Table  2,  the  deviations  of  the  calculated  from  the 
test  results  for  the  steady  rate  curves  are  generally  less  than  the 
experimental  accuracy  of  5%.  The  percentage  errors  in  comparing 
respective  stress  increment  values  are  greater,  but  the  overall 
characteristics  of  the  calculated  curves  for  sudden  strain  rate 
changes  are  similar  to  those  of  the  tests. 

Preference  in  matching  the  respective  steady  rate  curves  was 
given  to  the  room  temperature  case  which  influenced  the  constants 
in  the  Z (W  ) function.  The  calculated  steady  rate  at  523°K  also 
matched  fairly  well,  but  those  at  the  lower  temperatures  deviated 
fron  those  of  the  tests  by  more  than  5%  at  strains  higher  than 
about  15%.  For  this  reason,  only  results  to  15%  strain  are  shown 
for  77°K.  In  retrospect,  it  seems  that  somewhat  better  overall 
matching  could  have  been  achieved  with  a slightly  higher  value  for 
Zi  than  that  given  in  Table  1,  that  is,  Zi  should  have  been  about 
250  MPa.  However,  the  results  shown  match  the  test  curves  fairly 
well  are  are  considered  adequate  to  represent  the  potentialities 
of  the  analytical  representation.  i 
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Fig.  7 provides  a detailed  illustration  of  the  calculated 
changes  in  the  stress  and  the  strain  rate  components  consequent 
to  an  instantaneous  jump  of  six  decades  in  the  applied  strain  rate. 
The  example  shown  in  Fig.  7 is  for  the  second  jump  (•y^=12%)  at 
298°K.  Prior  to  the  jump,  at  the  low  strain  rate  Yg»  the  deforma- 
tion increments  are  almost  entirely  plastic  (99.2%).  Since  the 
plastic  strain  rate  y (=d'^)  is  a function  of  stress  which  doesn't 
change  instantaneously,  the  discontinuity  in  total  strain  rate  must 
be  accomodated  by  the  elastic  strain  rate  y which  becomes  essen- 
tially equal  to  the  imposed  rate.  This  results  in  an  elastic  stress 
increment  and  a rapid  increase  in  the  stress  level.  The  plastic 
strain  rate  correspondingly  increases  rapidly  and  becomes  almost 
equal  to  the  new  imposed  rate  in  about  lys  (Fig.  7) . During  this 
time  interval,  the  elastic  strain  rate  decreases  to  less  than  1% 
of  the  total.  At  all  the  test  temperatures,  the  calculated  re- 
adjustment in  stress  and  strain  rates  was  completed  in  less  than 
2ys  which  was  taken  as  the  jump  time  interval  for  obtaining  Ax. 

Since  the  actual  jump  times  in  the  experiments  [7]  were  about  lOys , 
the  initial  sharp  increment  in  stress  would  be  slightly  reduced  and 
would  be  followed  by  steady  plastic  flow.  Plastic  work  during  the 
strain  rate  jump  would  be  negligible  so  the  parameter  that  repre- 
sents the  microstructural  state,  Z,  would  remain  constant  during 
the  jump.  The  magnitude  of  the  stress  increment  could  then  be  cal- 


culated directly  from  (5)  for  the  appropriate  value  of  Z.  In  prac- 
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It  is  noted  that  if  the  strain  rate  jump  were  applied  locally  ; 

to  one  end  of  a bar,  then  an  elastic  wave  would  propagate  due  to  l 

the  elastic  deformation  increment.  This  is  consistent  with  the  ■ 

known  result  that  a stress  increment  suddenly  applied  to  a plasti- 
cally deforming  rod  would  propagate  at  the  elastic  bar  velocity. 

It  is  also  of  interest  to  follow  the  details  of  the  strain 
rate  jump  on  the  Df,  plane.  Fig.  1.  With  increasing  plastic 
deformation  at  the  lower  strain  rate,  the  hardness,  Z,  increases 

so  that  the  state  curve  shifts  on  the  D^,  plane  (point  A moving 
to  B corresponding  to  Z'  changing  to  Z").  At  B the  total  strain 
: 3 suddenly  increased  to  so  the  plastic  deformation  rate  rises 
in  a very  short  time,  about  2ms,  to  point  C corresponding  to  the 
increase  in  stress,  AJ^ [= ( A i j ‘ J , at  the  same  state  value  Z".  Con- 
tinued deformation  at  the  high  strain  rate  generates  additional 
plastic  flow,  C to  D,  v/ith  increasing  values  of  Z.  The  steady  high 
strain  rate  curve  would  have  state  points  corresponding  to  C and  D, 
namely  C and  D' , at  the  same  stress  and  plastic  strain  rates.  The 
extent  to  which  the  actual  incremental  curve  CD  matches  the  steady 
curve  C'D'  can  be  viewed  as  a measure  of  the  adequacy  of  a single 
internal  state  variable  in  the  constitutive  equations.  Examining 
the  details  of  the  stress-strain  curves  for  copper  in  [7]  on  this 
basis,  the  use  of  a single  microstructural  parameter  seems  reason- 
able but  is  not  precise.  However,  other  results  from  decremental 
temperature  tests  [11]  tend  to  support  the  applicability  of  a 
single  internal  state  variable.  It  is  noted  that  point  F in  Fig.  1, 
which  can  be  reached  at  a steady  high  rate,  has  the  same  strain  as 
points  B and  C but  corresponds  to  a different  state. 

- ^ 
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Another  measure  of  the  validity  of  the  representation  is  the 
comparison  of  the  magnitude  of  the  calculated  stress  increments  Ax 
with  those  obtained  in  the  experiments.  These  are  listed  in  Table  2. 
The  percentage  differences  are  greater  than  those  between  the  calcu- 
lated and  experimental  steady  rate  curves,  but  the  form  of  dependence 
of  the  calculated  stress  increments  with  strain  and  temperature  cor- 
responds to  those  of  the  tests.  As  in  the  tests,  the  calculated 
stress  increment  increases  with  increasing  strain,  7^.  The  tempera- 
ture dependence  of  the  Ax  values  is  erratic  at  the  lower  temperatures, 
but  the  calculated  results  do  follow  the  trends  of  the  tests.  Both 
sets  of  results  for  the  stress  increments  are  maximum  at  523°K. 

The  calculated  stress  increments  are  higher  than  those  of  the 

tests  by  15  to  70%  based  on  average  test  values.  There  are  a few 

reasons  for  these  differences  apart  from  using  the  steady  rate 

curves  as  the  basis  for  the  material  constants.  These  are:  the 

strain  rate  for  the  dynamic  incremental  tests  were  usually  lower 

• -1 

than  the  nominal  y^=300  sec  on  which  the  calculations  were  based; 
the  2ys  time  interval  used  for  the  calculated  Ax's  included  some 
work  hardening  plastic  flow;  the  experimental  time  for  the  change 
of  strain  rate,  about  lOys,  would  tend  to  diminish  the  initial 
stress  jump  which  was  taken  as  the  experimental  Ax . All  these 
factors  are  in  the  direction  of  increasing  the  calculated  stress 
increments  relative  to  the  observed  ones.  The  first  two  are  fairly 
small  but  the  third  could  be  of  some  significance.  These  factors 
might  account  for  most  of  the  differences  in  results  except  for  the 
largest  ones  at  523'’K.  Further  agreement  on  the  stress  increments 
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could  have  been  obtained  at  the  expense  of  poorer  matching  of  the 
steady  rate  curves.  Reducing  the  calculated  stress  increments 
is  achieved  mainly  by  increasing  the  n values,  i.e.  lowering  the 
rate  sensitivity.  A higher  n value  at  523°K  would  actually  permit 
better  fitting  of  a hyperbolic  relation  to  the  n(T)  function,  Fig. 
9,  i.e.,  eq.  (10b)  would  then  have  greater  validity. 

The  equations  for  (T)  and  n(T),  (9)  and  (10),  obtained  from 
the  results  of  the  numerical  exercises  follow  the  expected  trends 
that  hardness  would  decrease  and  rate  sensitivity  would  increase 
at  higher  temperatures.  That  a direct  correspondence  seems  to 
exist  between  plastic  working  and  lowering  the  temperature  is  an 
interesting  result  that  could  be  checked  experimentally. 

It  is  difficult  to  relate  eqs.  (9)  and  (10)  directly  to  para- 
meters in  the  rate  equations  for  plastic  flow  based  on  idealized 
dislocation  mechanisms.  The  exponential  form  (5)  could  be  approxi- 
mated by  a power  law  over  part  of  its  range  where  the  parameter  n 
would  be  equal  to  the  power  q multiplied  by  a constant.  Campbell 
shows  in  [10]  that  temperature  dependence  could  be  approximated 
in  such  a power  law  for  constant  rate  conditions  by  making  q=q'/kT. 
This  would  correspond  to  a hyperbolic  relation  for  n(T)  such  as 
indicated  by  (10b).  A strain-rate,  temperature  correspondence  in 
the  stress-strain  relation  is  also  implied  by  eqs.  (10),  but  this 
enters  in  a more  complicated  manner  than  the  simple  formula  pro- 
posed in  [9 ] . 

Results  given  in  [7,8,12]  for  the  metals  aluminum,  magnesium, 
zinc,  and  titanium  subjected  to  rapid  changes  in  strain  rate  show 
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characteristics  that  are  basically  similar  to  those  for  copper. 
Under  these  circumstances,  the  constitutive  equations  described 
in  this  paper  may  be  applicable.  For  certain  metals,  e.g.  mild 
steel,  molybdenum,  and  niobium,  the  strain  rate  jumps  sometimes 
lead  to  stress  levels  that  exceed  the  steady  rate  values  [12,13] . 
These  unusual  results  may  be  due  to  strain  ageing  effects  at  the 
lower  strain  rates  which  can  be  incorporated  into  the  constitutive 
equations . 
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List  of  Captions 


Fig.  1 - Plots  of  Da  = Fi(J2,Z)  for  various  hardened  states  Z 
showing  path  for  sudden  jump  in  strain  rate. 


Fig.  2 - Plot  of  hardness  parameter  Z for  copper  as  a function 
of  plastic  work  (per  unit  volume) . 

Fig.  3 - Calculated  and  experimental  static,  dynamic  and 

incremental  stress-strain  curves  for  copper  at  77°K. 

Fig.  4 - Calculated  and  experimental  static,  dynamic  and 

incremental  stress-strain  curves  for  copper  at  148°K. 

Fig.  5 - Calculated  and  experimental  static,  dynamic  and 

incremental  stress-strain  curves  for  copper  at  298°K. 

Fig.  6 - Calculated  and  experimental  static,  dynamic  and 

incremental  stress-strain  curves  for  copper  at  523°K. 


Fig.  7 - Calculated  time  variation  of  stress  and  strain  rate 
components  due  to  instantaneous  change  of  applied 
strain  rate  from  2 x 10  “^sec  ^ to  3 x lO^sec 


Fig.  8 - Dependence  of  initial  value  of  hardness  parameter,  Z^, 
on  test  temperature. 


Fig.  9 - 


Dependence  of  strain  rate  sensitivity  parameter,  n, 
on  test  temperature. 
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TABLE  1 


1. 


Constants  in 

Z = Zi 

m = mo 


Z (Wp)  function  for  copper 

+ ( Zo -Z : ) exp ( -mW  ) 

P 

+ mi  exp ( -aWp) 


(temperature  independent) : 

(6) 

(7) 


Zo  = 31  (MPa) 

Zi  = 237  (MPa) 
mo  = 1500  (MPa) 
m,  = 2500  (MPa) 
a = 5000  (MPa) 


2,  Temperature  dependent  material  constants: 


T°K 

Z^(MPa) 

I 

n 

77 

82 

14. 3 

148  j 

49 

13.4 

298 

31  , 

, 9.2 

1 

523 

1 26 

3.6 

w I 

( 

I 

! 3.  Temperature  independent  material  constant: 


j 

1 


I 

I 


TABLE  2 - COMPARISON  OF  RESULTS 


Note:  bars  indicate  average  value  of  several  tests;  Y-  = strain  at  jump;  other  symbols  defined  in  Fig. 


(DdW)  Z 


(MPa) 
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— experimental  [7] 
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T=  2980K 

O^  = 3x10^  sec’ 

% = 2x10'^  sec"’ 

experimental  [?] 
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